On the Stability of Circular Orbits of Particles Moving around Black Holes 
Surrounded by Axially Symmetric Structures 
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The Rayleigh criterion is used to study the stability of circular orbits of particles moving around 
static black holes surrounded by different axially symmetric structures with reflection symmetry, 
like disks, rings and halos. We consider three models of disks one of infinite extension and two finite, 
and one model of rings. The halos are represented by external quadrupole moments (either oblate 
or prolate). Internal quadrupole perturbation (oblate and prolate) are also considered. For this 
class of disks the counter-rotation hypothesis implies that the stability of the disks is equivalent to 
the stability of test particles. The stability of Newtonian systems is also considered and compared 
with the equivalent relativistic situation. We find that the general relativistic dynamics favors the 
formation of rings. 

PACS numbers: 04.70.Bw, 95.10.Eg, 98.35.Mp, 98.62.Mw 



I. INTRODUCTION 

One of the main tools used to study the gravi- 
tational field of any distribution of mass, e.g., a black 
hole (BH) , is to examine the motion of nearby test parti- 
cles. In other words, the complete understanding of the 
different solutions of the geodesic equation leads us to 
the comprehension of the gravitational properties of the 
center of attraction, in the present case a black hole sur- 
rounded by and axially symmetric structure. Because of 
its mathematical simplicity and physical relevance circu- 
lar orbits play a particularly important role in the study 
of the gravitational properties of a central body. Also 
circular orbits usually represent limit situations between 
different regimes. 

The purpose of this paper is to study the stability 
of circular orbits in different gravitating systems formed 
by a BH surrounded by axially symmetric structures. In 
particular, we study orbits in BH + disks systems, BH + 
rings systems, and BH +multipolar fields. In the case of 
orbits inside the matter, e.g. a disk, we have that in gen- 
eral the stability of the circular orbits can be associated 
to the stability of the matter itself. The method used 
to study the stability is the Rayleigh criteria of stability 
(RCS) for a rotating fluid that can be easily adapted 
to the study of orbits of test particles. In the context 
of general relativity the RCS has been used to study the 
stability of different classes of self-gravitating disks, see 
for instance 

The paper is divided as follows: In Sec. II we 
review the RCS and in Sec. Ill the Einstein equations 
for the superposition of a static BH and an axially sym- 
metric structure with reflection symmetry is presented. 
In Sec. IV the stability of circular orbits moving around 
the superposition of a BH and three models of disks are 
considered, one of infinite extension and two finite. In 
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Sec. V we study the stability of particles orbiting on the 
plane of a ring + BH system, inside as well as outside 
the ring. In Sec. VI we consider circular orbits for a 
BH perturbed by either an internal quadrupole field or 
a external quadrupole field. We consider oblate as well 
as prolate quadruples. We have that some general fea- 
tures found in the study of the previous examples can be 
simulated with these quadrupole fields. In the last sec- 
tion we discuss some of the previously found results. We 
also review in an Appendix the so called counter-rotating 
hypothesis. 



II. RAYLEIGH CRITERIA OF STABILITY 

Consider, in Newtonian dynamics, a small body 
of mass \i moving in a circular orbit of radius r = ro 
around a fix center of gravitational attraction. In the 
non inertial frame associated to the moving body we have 
an equilibrium situation: the gravitational force [fg^o)] 
is equal to the centrifugal force [F c (ro) — L 2 (r )/ (fir^)], 
where L is the angular momentum of the particle of mass 
(i. Now let us virtually displace the moving particle to 
a higher orbit (r > ro) keeping the same angular mo- 
mentum. Then the centrifugal force in this new posi- 
tion is F c (r) — L 2 (ro)/(/ir 3 )]. To have a situation of 
equilibrium the displaced particle should move down- 
ward, i.e., in direction to the initial orbit. In other 
words, the gravitational force in r should be greater 
than the centrifugal force, F c (r), i.e., F g (r) > F c (r), but 
F g (r) = F c (r) = L 2 (r)/(/w 3 ). Hence L 2 (r) > L 2 (r ), 
by doing a Taylor expansion of L 2 (r) around r = ro we 
find that, > 0, for a stable circular orbit. Rayleigh 
in the original derivation of his criterion [JJ , considers a 
rotating ring of fluid under the force of a pressure gra- 
dient instead of a the particle moving in a circular orbit 
attracted by a gravitational field. 

Now let us consider a particle moving in a static 
background metric with axial symmetry in spherical co- 
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ordinates, (t,R, i?, <p), 

ds 2 = g tt dt 2 + g RR dR 2 + gmd-d 2 + g w dip 2 , (1) 

where the metric functions g^ v depends only on the vari- 
ables i?, ■&. The geodesic equation for a circular motion 
on the plane z? = 7r/2, gives us the motion equation, 



gtt,ni 2 + g vv ,B.<f> 2 = 0, 



and the constants of motion, 



1 



gut 2 +g^ 2 , 



E = gut, 

h = g<p<p<p, 



(2) 

(3) 
(4) 
(5) 



where the orverdot denotes derivative with respect to the 
proper time s and gu,R = dugu, etc. The constant E rep- 
resents the relativistic specific energy and h the specific 
angular momentum. 

Note that the motion equation @ can be cast as 
a balance equation valid on the plane d = 7r/2 , 



9u,rE 2 jg\ t = -g vv ,nh 2 / 'g 2 v . 



(6) 



So as in the Newtonian case we have a balance between 
the "gravitational force" and the "centrifugal force". 
Thus, assuming the metric represents the gravitation of 
a central body we will have stability of circular orbits on 
the plane z? = tt/2, when 



hh r > 0. 



From (J2J- @we find, 
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g vv gu,R 



gug vv ,R - gu.Mw 



E 2 



g u g<pv,R 



gttg V ip,R — gu,Rg v <p 



(7) 

(8) 
(9) 



Since gravity is an attractive force we can repeat 
the same reasoning that leads to 10, but now with the 
constant E(R, tt/2) instead of h(R, tt/2), we find 



EE, R > 0. 
From © and © we get, 



hh 



9tt 



EE 



,R- 



(10) 



(11) 



But, always —g vv /gtt > 0, then the relation l(TU|l gives us 
the same information than (Q, in other words we have 
consistence. 

For the stability of circular orbits in stationary 
axisymmetric spacetimes see Ref. [j| and for stability of 
rotating fluids in a relativistic stars see Ref. . This last 
reference presents a very general study of the stability 
of relativistic fluids that generalizes the Lord Rayleigh 
work by considering a relativistic fluid with both bulk 
and dynamical viscosity and also heath flow. A slightly 
different treatment of the stability of circular orbits like 
the one presented here can be found in Ref. 7]. 



III. WEYL METRICS: BLACK HOLES + 
AXIALLY SYMMETRIC STRUCTURES 

The external gravitational field produced by an 
axially symmetric body can be well described by the Weyl 
metric p|, 

ds 2 = e 2 ^dt 2 - e- 2 ^[r 2 d V 2 + e 2 ^(dr 2 + dz 2 )}, (12) 

where the functions ijj and 7 depend only on the cylindri- 
cal coordinates r and z. The vacuum Einstein equations 
[Rau — 0) reduce to the usual Laplace equation in cylin- 
drical coordinates, 

1pj-r + 4>,r/r + Ip.zz = 0, (13) 

and the quadrature, 

tyty] = r[(ip 2 r - ^ 2 z )dr + 2ip tr il> tZ dz}. (14) 

When ip satisfies the Laplace equation this differential is 
exact. 

The Schwarzschild solution in Weyl coordinates 
takes the form, 



1, R+ + i?_ - 2m 

W Q = — 111 , 

v 2 R++R_+2m 

1, (R++R) 2 -4m 2 
Is = — hi . 

' 2 R+R- 



where 



R ± = yjr 2 + (z ± to) 2 



(15) 
(16) 

(17) 



The function ips is just the Newtonian potential of a bar 
of length 2m and density 1/2. The Weyl coordinates has 
the drawback that the Schwarzschild metric does not look 
spherically symmetric and that the horizon is squeezed 
into a line of length 2m. But has the advantage that one 
of the Einstein equations is the Laplace equation that 
is linear, hence it allows to consider different composite 
systems like a black hole surrounded by disks or halos in 
an exact way. 

The coordinate transformations, 

r= [R(R - 2m))]i sm<&, z = {R — m) cos (18) 

takes the black hole metric with (|15[l - Hltj|) m t° the 
usual Schwarzschild form. 

We shall consider solutions of the form, 



-0 = ips + i>, 



(19) 



where ip{r, z) is a solution of either Laplace equation or 
Poisson equation with a density that has support on the 
plane z — (t? = 7r/2). The function 7 in this case is 
also written as 



ibP] =ls + l(r, z). 



(20) 
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In the case of vacuum solutions j(r, z) can be computed 
directly from {Tilt and [t*H3|) . 

Therefore for the superposition (|19[) the Weyl met- 
ric takes the form, 



and its derivative can be cast us, 



ds 2 = (1 _ ^yU dt 2 
R 



e 2(7-0) 

-i 2 m 

1 H 



di? 2 



-R 2 e~ 2 " i '(e 2Al M 2 + sin 2 r3aV> 2 



(21) 



When ip = 7 = we recover the Schwarzschild metric. 

We shall analyze three types on superpositions of 
BH and matter. In the first, the matter is represented by 
a potential ip that describes a thin disk. In this case we 
will have that only two components of energy-momentum 
tensor (EMT), T\ and Tg, are different from zero on the 
plane of the disk. Thus we have a disk with no radial 
pressure or tension to support the gravitational attrac- 
tion. Since the Weyl metric is static we have no rotation. 
Then a counter-rotation hypothesis is needed to have a 
stable situation: we assume that on the disk as many par- 
ticles move clockwise as counterclockwise (see Appendix 
A). Even though, this interpretation can be seen as a 
device, there is observational evidence of disks made of 
streams of rotating and counter-rotating matter 9]. We 
believe that counter rotating hypothesis appeared first, 
in this context, in Ref. fl(ij ]. 

In the general case we will have radial pressure 
and magnetic forces and our simpler analysis need to be 
modified. In the case when no magnetic fields are present 
and the density and pressure of the matter of the disk are 
low, e.g., a diluted gas of stars, the main contribution to 
the force on each fluid particle is the gravitational force 
that is taken into account through the spacetime metric. 

In the second case we study the superposition of a 
BH and a thin ring. We consider orbits inside the ring as 
well as orbits outside the ring. And in the third case we 
consider the superposition of a BH with multipolar fields. 
We first consider internal multipolar perturbations, i.e., 
if) is a multipolar field that is zero at spatial infinite. In 
this way we can describe the gravity outside compact 
systems formed by a true black hole (or a dense object) 
surrounded by a distribution of matter like a ring or a 
small disk. Also we can have an axially symmetric static 
dense object with either polar deformations or polar jets. 
We also consider the case of a BH with external multipo- 
lar perturbations, i.e., ip is a multipolar field, solution of 
Laplace equation (|13|l . that diverges at spatial infinite. 
In this case we can represent the gravity in the region 
limited by a black hole and an outer axially symmetric 
shell of matter or halo. 

All the models of BH with surrounded matter that 
we consider in this paper are constructed from superpo- 
sition of the Weyl solutions of the form H2U|) and its asso- 
ciated spacetime is described by the metric (|21|) . In this 
case the constant of motion JjJ takes the form, 



h 2 = 



R 2 e- 2 ^[m + R(R - 2m)f fl ] 
R-3m-2(R-2m)i> tR 



(22) 



hh o = 



Re- 2 v[m(R- 6m) +F 1 ] 

2(^- 3771 - 2^2) 



where 



Fx = i?(28m 2 - ISmR + 3i? 2 )^ jfl 
6R 2 (R - 3m) (R - 2m)^ 2 R + 
+R 2 (R-m){R-2m)^ :RR , 

F 2 = R(R-2m)i) tR 



4i? 3 (i? 



(23) 



2m) 2 ^ R 



(24) 
(25) 



In the present case the motion is limited to the plane z — 
0(r3 = tt/2). From fTJfy we have that r = [R(R - 2m))]i 
on the plane and t£(r,0) = 4>([R(R - 2m))] 3,0). Note 
that when tp = 0, pure BH case, from ill'MI) and hh^n = 
we have R = 6m, as required. This orbit is known as 
the last stable circular orbit (LSCO) or the marginally 
stable circular orbit (MSCO). In our study of stability we 
will use Schwarzschild like coordinates. The use of Weyl 
coordinates or proper distance coordinates are known to 
give similar qualitative results [Tl| . 



IV. BH + DISKS 

In this section we study the superposition of a BH 
and three different models of disks. In the first case the 
disk is the relativistic generalization of the Newtonian 
Kuzmin disk that is studied in Ref. |2|. This solution 
is closely related to the vacuum solution of the Einstein 
equation known as Chazy-Curzon metric |12|. The metric 
function ip for the superposition of this disk with a BH 
is, 



$(r,0) = -M/(r 2 + a 2 )3 



(26) 



The disk has mass M an infinite radius, a is a parameter 
related to the maximum of the energy density. 

The second and third model of disks are the zeroth 
and first Morgan and Morgan disks [ljj that are obtained 
by solving the Laplace equations in oblate coordinates. 
The function ip(r,Q) for the zeroth order Morgan and 
Morgan disk is 



- M a 

tp{r, 0) = arctan -===. 

a \/r 2 - a 2 



$(r,0) 



nM 
~2a 



■, r < a. 



r>a (27) 
(28) 



M represents the mass of the disk and a its radius. This 
disk has infinite density on the edge. The superposition 
of this disk and a BH is studied in [l3|. 

For the first Morgan and Morgan disk we have, 

= — ^— ^[V r 2 — a 2 — a _1 (r 2 — 2a 2 ) arctan —== 



4a 2 
3?rM . 



(a 2 -7-72), r<a 



\Jr 2 — a 2 
r> a. (29) 

(30) 



4 



Again M is the mass of the disk and a its radius. The su- 
perposition of this disk with a BH was studied in 0] . 
This disk has finite density on the rim. For simplic- 
ity we have presented the potentials Weyl coordinates 
(t,r,z,(p). We note that to a disk of radius a in Weyl 
coordinates corresponds a disk of radius 



R a = m+ ^Jm? + a 2 , (31) 

in the usual BH coordinates that will be used in our anal- 
ysis of stability. We shall study numerically the condition 
hhn = that will give us information on the stability of 
circular orbits for the different superpositions of BH and 
disks. 

For the Kuzmin-Chazy-Curzon disk in Fig ^ we 
study the curve hh_R — for different values of the pa- 
rameter a, we take a = 0.5, 2.5, 4.5, 6. Along all this work 
we take units such that G = c = 1 and we fix the unit of 
length by taking the mass of the central BH as one. We 
have stability in the regions limited by the upper curves 
and the axis R and in regions limited by the lower curves 
and R = 2m. In the first case, a = 0.5, for a given total 
mass parameter of the disk, M, we have a region, close to 
the BH horizon, wherein we have stable circular orbits. 
For larger values of R we have a region of instability. In 
principle this region of instability will give rise to a gap in 
the disk. For even larger values of R we will have another 
region of stability that extent to infinite. In the second 
case, a — 2.5, closer to the horizon, for 4 < M < 6, will 
appear another region of instability. So we can have the 
formation of stable flat rings in this region. In the third 
case, a = 4.5 we can have rings around M = 3 and for 
M > 8. In the last case, a = 6, rings can also be formed. 
It is interesting to compare this case with the equivalent 
Newtonian situation. By applying the RCS to the system 
represented by the Newtonian potential 

<f> = -M/(R 2 + a 2 )2 - m /R (32) 

we find only stable circular orbits. Therefore the forma- 
tions of rings and gaps, in the present case, is a pure 
general relativistic effect. 

Now we shall consider the zeroth order Morgan 
and Morgan disk. In the interior of the disk we have 
that the function -0 is constant. Hence the condition 
hhn = gives us R = 6m, the LSCO of the usual BH. 
Therefore if the radius of the disk is less than 6m it is 
completely unstable. For a disk of radius grater than 6m 
we will have that the region of the disk between 6m and 
its radius will be stable. The condition hh j) — outside 
the disk is studied in Fig|2| 

We have stability in the region limited by the up- 
per curve and the axis R and in the region limited by 
the lower curve and the horizontal line that represents 
the radius of the disk. Therefore we can have the forma- 
tion of a stable structure orbiting around the disk close 
to its edge. After this zone of stability we find a zone of 
instability followed by one of stability. 



BH + KCC (a=0.5) 




1 2 3 4 5 6 

M 

BH + KCC (a=0.5) 




2 4 6 8 10 

M 

BH + KCC (a=4.5) 




5 10 15 20 

M 

BH + KCC (a=6) 




I 1 1 1 1 

5 10 15 20 

M 



FIG. 1: Curves hh t u — for the superposition of a BH 
and a Kuzmin-Chazy-Curzon disk for different values of the 
parameter a. The curve, R — 2m, is also shown. In the first 
three cases we have stability in the regions limited by the 
upper curves and the axis R and in the regions limited by the 
lower curves and R = 2m. In the last case we have stability 
in the region limited by the upper curve and the lower curve. 
We use units such that G — c — m — 1. 
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BH + MMO (a=10) 



BH + MM1 ext. (a=9 ,12) 





FIG. 2: Curves hh t R — for the superposition of a BH and 
a zeroth order Morgan and Morgan disk for a — 10. The 
horizontal line indicates the position of the disk boundary. 
We have orbit stability in the region limited by the upper 
curve and the axis R and in the region limited by the lower 
curve and radius of the disk. The units are the same as in the 
previous figure. 

BH + MM1 (a=6,9,12) 




FIG. 3: Curves hh t R — for the superposition of a BH and 
the first Morgan and Morgan disk. The top curve represent 
a disk with parameter a — 12, the middle with a = 9 and the 
bottom curve with a = 6 We have stable circular orbits in the 
region above each respective curve. 



It is instructive to compare these results with the 
equivalent Newtonian situation that is represented by the 
potential, 



M a 
arctan — ■ 



rn 
~R' 



(33) 



We find a very different situation, we only have a small 
zone of instability close to the disk boundary. The inte- 
rior of the disk is stable. 

Now we shall consider the superposition of first 
Morgan and Morgan disk with a BH. In Fig^Jwe present 
the curve hh t R = for three different disks with parame- 
ters a = 12 (top curve), a = 9, and a = 6 (bottom curve). 
We have a region of instability in the region under the 
curve hhj-i = in each case. We have that the presence 
of a disk lower the black hole's LSCO. 

The curves hh^ = for the superposition of a 
BH and the first Morgan and Morgan disk with a = 9 



FIG. 4: Curves hh,R = for the superposition of a BH and 
the first Morgan and Morgan disk with a = 9 (bottom curve) 
and a — 12 (top curve). We have stable circular orbits in the 
region limited by the curve and the axis R. 



is presented in Fig. |2]for a region outside the disk. We 
have that the closer to the edge of the disk we always 
have unstable circular orbits. For larger radius we have 
stable orbits. 

It is interesting to compare this case with the same 
Newtonian situation. For the interior of the disk with 
a Newtonian center of attraction we find no regions of 
instability. Outside the disk the situation is similar to 
the one shown in Fig. but in the Newtonian case the 
curves are lower. 

The superposition of the BH with the zeroth order 
Morgan and Morgan disk and with the first Morgan and 
Morgan disk present different properties of stability near 
the edge of the disk. We have that the first mentioned 
disk has a singular rim and the second has a regular edge. 
This last case represent a physically acceptable situation. 



V. BH + RINGS 

In this section we consider the superposition of a 
black hole and a family of thin rings. The function %p for 
the thin rings will be taken as, 



1 



2 Wr 



arctan ■ 



M , 



2 l Va? 



r < a. 



> <34) 
(35) 



The parameters M and a represent the mass and radius 
of the ring, respectively. This family of thin rings as well 
as its superposition with a black hole is studied in some 
detail in [l]|. 

In Fig. we consider the curve hh^R = for the 
superposition of a BH and the ring with potential <j> given 
by l|34() - (|35|l . The first graph shows the region exterior 
to the ring and the second the interior region for a = 8 
(left hand side curve) and a = 10 (right hand side curve). 
For a = 8(a = 10) we have R 8 - 9(iZi - 11). 
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BH + Ring (a=8) 




2 4 6 8 10 

M 

BH + Ring (a=10) 




10 1 ' ' ' ' 1 

2 4 6 8 10 



M 

BH + Ring int. (a=8,10) 




0.25 



M 



FIG. 5: Curves hh,R — for the superposition of a BH and 
a thin ring. The first two graphs show the region exterior to 
rings with a = 8 and a = 10 and the third the interior region 
for rings of the same radii the LHS (RHS) curve correspond 
to a = 8(10). In the first two graphs we have stable circular 
orbits in the the regions limited by the upper curves and the 
axis R and in the regions limited by the lower curves and the 
horizontal line that represents the location of the ring. In the 
third graph we have stable circular orbits in the region limited 
by the corresponding curve and the axis R. 



For the exterior region we have stable circular or- 
bits for values of (M, R) limited by the upper curve and 
the axis R and in the region limited by the lower curves 
and the horizontal line that represents the location of the 
ring. In the third graph we have stable circular orbits in 
the region limited by the corresponding curve and the 
axis R. 

The regions of stability for the equivalent New- 



tonian situation are different. We have only one small 
region of instability of circular orbits with radii larger to 
the ring radius. For the interior of the ring we have that 
the curve hh y n — is similar to the precedent case for 
values of R close to the ring radius. The Newtonian case 
is similar to the relativistic case, but the curve does not 
turns back to R = 6m it decreases as R ~ 1/M 1 / 3 . The 
analysis of stability of circular orbits for the superposi- 
tion of a BH and a flat ring (Lemos-Letelier solution |14|h 
as well as other aspects of this solution were considered 
in great detail in [lij, [ill, an d E3- 

VI. BH + MULIPOLAR FIELDS 

Multipolar exact perturbation of a BH can ap- 
proximate all situations studied in the previous sections 
of this article that refer to compact sources of gravity. 
Also we can have some insight for the stability of circu- 
lar orbits in the generic case of a static BH with axially 
symmetric perturbations with reflection symmetry. 

The functions tp( r : z ) i n the present case are 

i>(r,0) = -QJr 3 (36) 
?/>(r, 0) = -Q e r 2 (37) 

Where Qi represents the quadrupolar moment of an in- 
ternal distribution of matter, i.e., of sources of the gravi- 
tational field near the BH, e.g., disks, rings, etc. This 
field will approximate the gravitation field outside its 
source. We have two cases Qi > and Qi < 0. The 
first case represents oblate deformations like disks and 
rings and the second prolate deformations, e.g., jets. The 
constant Q e is related to the quadrupolar moment of ex- 
ternal distributions of matter, i.e., of sources of the grav- 
itational field far away from the BH, e.g, the gravity in 
the space limited by a ring or a shell of matter (halo). 
We also have two different situations in this case Q e > 
and Q e < 0. The first case describes oblate deforma- 
tions like a ring or oblate halos and the second prolate 
deformations like prolate halos. 

In Fig. we present curves hh t R = for the 
superposition of a BH and an internal multipolar field 
with <j) give by (|36|l . The first graph shows the oblate case 
and the second prolate deformations. In the oblate case 
we have stable circular orbits in the the regions limited 
by the upper curve and the axis R and in the regions 
limited by the lower curves and the horizontal line that 
represents the BH horizon. The most relevant region 
is the one with R > 6m. We see that the region of 
stability increases with the internal quadrupole moment. 
It is interesting to compare this result with the previous 
case. Let us consider the superposition of a BH with 
the first Morgan and Morgan disk. We see in Fig Q] that 
for a given disk mass the region of instability increases 
with the disk radius the same is true for rings, see Fig 
We also checked this property for the zeroth order 
Morgan and Morgan disk. We have, for these compact 
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BH + Quadrupole (int.) 
12 i 1 1 1 




I 1 1 1 1 1 

1 2 3 4 5 



Qi 

BH + Quadrupole (int) 




-1 -0.9-0.8-0.7-0.6-0.5-0.4-0.3-0.2-0.1 



Qi 

FIG. 6: Curves hh^n = for the superposition of a BH and an 
internal quadrupolar field. The first graph shows the oblate 
case and the second the prolate deformations. In the first we 
have stable circular orbits in the the regions limited by the 
upper curve and the axis R and in the regions limited by the 
lower curves and the horizontal line that represents the BH 
horizon. In the second graph we have stable circular orbits in 
the region limited by the axis R and the curves. 

objects, that the quadrupolar moments scale as Qi ~ 
Ma 2 . Therefore to a bigger radius corresponds a bigger 
quadrupole moment and a larger region of instability. 

In the prolate case, last graph of Fig. [5] we see 
that the quadrupole moment increases the zone of stabil- 
ity, with no quadrupole moment we have stability only 
when R > 6m. This is a rather surprising effect valid 
only for equatorial circular orbits. For generic orbits, a 
prolate internal deformation introduces great instability, 
moreover we can have very irregular orbits (chaos) |18| . 

In the equivalent Newtonian system, for the oblate 
case we have a region of instability near the center of at- 
traction. The limit curve is R = {\/SQi/m. We have no 
instability in the prolate case. Again for non equatorial 
circular orbits we can have chaotic motion in this case 




In Fig. [3 we present curves hh r = for the su- 
perposition of a BH and external multipolar field with <j) 
give by J37J. In the first graph, oblate case, we have sta- 
ble circular orbits in the the regions limited by the curve 
and the axis R. We have that this external quadrupole 
fields severely restrict the possibility of stable orbits. We 
have that for larger value of Q e > 0.00015 we do not 



BH + Quadrupole (ext.) 




2 I 1 1 1 1 1 

-0.1 -0.08 -0.06 -0.04 -0.02 
Qe 



FIG. 7: Curves hh,R = for the superposition of a BH and 
an external quadrupolar field. In the first graph , oblate case, 
we have stable circular orbits in the the regions limited by the 
curve and the axis R. In the second graph, prolate case, we 
have stable circular orbits in the region above the curve. 

have equatorial stable circular orbits. For Q e < 0.00015 
we have that the region of stability is limited from above 
as well as from below. For generic orbits in BH exterior 
halo systems we have chaotic orbits [P9| . 

We have that the exterior quadrupole moment of 
a ring or any external configuration of matter like a shell 
or halo scales as Q e = M/a 3 , with M the mass of the 
shell and a the typical size. Then large radius indicates 
small quadrupole moment and large stability region. We 
note that second graph of Fig. [S] illustrate exactly this 
conclusion. In the equivalent Newtonian case we have 
that the curve hh, R = is R — (m/^Qe) 1 ^ . 

In the prolate case, second graph of Fig. [7J we 
have stable circular orbits in the region above the curve. 
We see that the prolate external matter increases the 
region of stability. For the Newtonian equivalent system 
we always have stable orbits. Marginally stable orbits 
in quadrupole deformed centers of attraction were also 
studied in [2fj|. 

VII. CONCLUDING REMARKS 

We think that one of the most important results 
shown in this paper is that we can have the formations of 
rings and gaps around a BH as a pure general relativistic 
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effect, of course this is a consequence of the presence of 
the BH horizon that modifies the usual Newtonian effec- 
tive potential adding new critical points. These critical 
points change completely the dynamics of the orbiting 
particles. In particular, the maximum has the property 
of being a saddle point in phase space. Hence perturbed 
orbits tend to have irregular motions. 

The use of multipoles to represent a distribution of 
matter in the general relativistic context need to be taken 
with care. In general relativity, the multipolar expansion 
take into account that we have a tensorial source of the 
gravitational field and not a scalar one as in Newtonian 
gravity. We can have different sets of general relativistic 
multipoles that in the limit of Newtonian gravity give the 
same set of usual Newtonian multipoles, for a discussion 
of this point in the context of static axially symmetric 
spacetimes see |2l|. We found for orbits that are a few 
Schwarzschild radius away from the center of attraction 
no significant discrepancy due the use of different mul- 
tipolar expansion with the same Newtonian limit |2l| . 
Hence there is no great error in taken any multipolar ex- 
pansion with the same Newtonian limit, the one taken in 
this paper is the simpler. 

Our study of chaos in BH with halos an d de- 
formed centers of attractions indicates that the sta- 
bility of circular orbits is quite different from the stability 
of generic orbits. This is not a surprising result since the 
equatorial circular orbits are a very special case of highly 
symmetric orbits. 
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APPENDIX A: THE COUNTER-ROTATING 
HYPOTHESIS 

The EMT for the counter-rotating matter is the 
sum of the energy momentum tensor of two pressureless 



streams of particles, 

rpflU _ rp(J,U _|_ rpfJ.U (Al) 

where 

T£ v = p+ U% Ul , Tt v = p_ UH U1 . ( A2) 

p+ and p- are the matter density of each stream that are 
considered equal (p + = />_). U+ and Ut are the corre- 
sponding normalized streams' four-velocities (UttU+n = 
U^U-n = 1). These four velocities can be written as 

U» = aU t *+l3(t> v : U» = aU^ - p<t? , (A3) 

where L™ is a timelike vector and <f> v an spacelike vec- 
tor tangent to one of the two streams of particles, also 
U^U„ = -cj) u (i) v = 1, U»<j) v = 0. Note that a 2 - [3 2 = 1. 
Therefore the EMT ijHJ can we put in the form 

= pU^U v + p^<t> v , (A4) 

with p — 2a 2 p + and = 2(1 — a 2 )p. 

Thus the counter-rotating hypothesis is consistent 
with the disks represented by the Weyl metrics solutions 
to the Einstein equations that have only azimuthal pres- 
sure. Furthermore we can assume that 

V^Tf = V M T^ = 0. (A5) 

From we have U^V ^ = UtV \JJ V _ = 0, i.e., each 
stream follows a geodesic motion. Note that (|A5() implies 
that V^T' 1 " = 0. Thus the supposition that each stream 
follows a geodesic flow is consistent with the Bianchi iden- 
tity for the Weyl metric. For multifluids, see for instance 
Ref. 

In summary, when the counter-rotating hypothe- 
sis is assumed, the stability of circular geodesies is equiv- 
alent to the stability of the flow of each fluid component. 
For the counter-rotating model of cold disks we have that 
the LSCO defines exactly the inner radius of the disk. 
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